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The SU {N) non-linear a-model [Q dynamics can be described by using a 'poly- 
nomial' action which is first order in the derivatives of the fields |0, ^ ^. The 
basic idea is that the usual action may also be written in terms of a {SU{N)) 1- 
form field which satisfies the Maureer-Cartan (flatness) equation F^y{L) = ||^. 
The polynomial formulation consists in writing the original action in terms of the 
vector field, adding a Lagrange multiplier term which enforces the Maureer-Cartan 
equation. As F^^, is a 2-form field (i.e., an antisymmetric tensor), so must be the 
Lagrange multiplier. 

From the usual description based on fields with values on a non-linear group 
manifold, and with the proper group-invariant path integral measure, one goes to 
another one in terms of a non-Abelian vector gauge field plus an antisymmetric (2- 
form) tensor field, both with linear path integral measures In three and higher 
spacetime dimensions there appears a gauge symmetry under transformations of 
the antisymmetric tensor field ^ which, when gauge fixed, requires of course the 
introduction of the corresponding Faddeev-Popov ghosts. In four dimensions, this 
symmetry is reducible 0, thus a proper BRST treatment prescribes the introduction 
of ghosts of ghosts. Note also that the bosonic part of the action is first order in the 
derivatives, and, depending on the gauge fixing one uses, the ghosts' action may or 
may not be first order. 

The 0(2) (~ U{1)) non-linear a-model in two (Euclidean) dimensions (or X — Y 
model) is the simplest possible example of application for the polynomial formula- 
tion. Besides the non-existence of the gauge symmetry linked to the antisymmetric 
tensor field in higher dimensions, this field can be decomposed as: O^y = e^u 
where 6' is a pseudoscalar field. Despite these simplifications, the model itself is 
far from trivial, and in its Statistical Mechanics version it undergoes the celebrated 
Kosterlitz-Thouless phase transition P]. Indeed, this phase transition occurs pre- 
cisely because the system has room for the existence of vortices. It is also well known 
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that it can be mapped to a Sine-Gordon model P], the correspondence between the 
0(2) spin field and the Sine-Gordon field being non-local. Thus the question presents 
itself how does one introduce vortices or any other singular configuration. 

The object of this paper is to study the polynomial formulation for the 0(2) 
model in two dimensions, showing how some familiar results of the usual formulation 
reemerge in a simpler way, and also how some extensions can be implemented. In 
particular, the Sine-Gordon description is obtained simply by integrating out the 
vector field and thus deriving an 'effective' action for the Lagrange multiplier, which 
becomes precisely equal to the Sine-Gordon one. 

The structure of this paper is as follows: In section 1 we briefly review the 
model in its usual formulation, describing some properties which we reformulate in 
their polynomial form in section 2. In section 3 we construct the mapping to the 
Sine-Gordon model, and in section 4 we show how to introduce strings of vortices. 

1 The non-polynomial or 'second order' formula- 
tion 

The model is usually deflned in terms of the Euclidean action (we follow the pre- 
sentation of ref. 0) 

^ = ^ /c^'x d,s{x) ■ d,s{x) (1) 
where s(x) is a two-component, real, continuum spin fleld 

s{x) = (si(x),S2(x)) , s^(x) = 1 , (2) 

and t is a (dimensionless) parameter which plays the role of a temperature. 
To solve the constraint on the modulus of s one can parametrize it as 

s(x) = (cos 0(x), sin 0(x)) (3) 
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where now (j) must be pseudoscalar under parity transformations of the two- 
dimensional spacetime (reflections about one of the axis), and one should note that 
physical quantities are 27r-periodic functions of (i.e., it is an angular variable). 
The 0(2) symmetry has been transformed into invariance under rigid translations 
of 0. The partition function is then 

Z = |p0exp[-5(0)], (4) 

where 

^(0) = I d'x ^^d,<P{x)d,<P{x) . (5) 

Every correlation function of spin variables can be obtained by linear combination 
of correlation functions of exponentials of the field (p 

n 

{ n exp[iej(p{xj)] ) , ej = ±1 , Wj (6) 

where the the arguments of the spin fields in the corresponding correlation 

function. In the spin-wave approximation one neglects singular configurations re- 
lated to the periodicity of 0, and then the correlation function (6) can be exactly 
calculated as 

( n exp[ze,0(x,)] ) = (-) ' n I - \)^'^'' , (7) 
j=i ^ j<k 

where A and m are UV and IR cutoffs respectively. When one takes the m — >■ limit, 
the only non-zero correlation functions are the ones which satisfy the 'neutrality' 
condition: 

n 

E^. = 0' (8) 

and this is precisely the condition for 0(2)-invariance of the correlation functions; 
i.e., invariance under translations of 0. In particular this implies the vanishing of the 
average of the one-spin function (as prescribed by the Mermin- Wagner theorem |1I0| ) 
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and, regarding the 2-point correlation function, (|^) yields 

A 27r t_ 

(s(x) ■ s(?/)) = (— ) {m\x-y\) 2-. (9) 
m 

So far the discussion has been confined to the low-temperature phase, characterized 
by an algebraic decaying of the 2-point correlation function (|^). Let us consider 
now the high-temperature phase, where vortices can appear unconfined, and the 
periodicity of is important. The vortices are finite-energy configurations such that 
the angular field changes by 2?™, where n is an integer, when one moves around 
a point, the 'center' of the vortex, n is known as the 'winding' of the vortex. A 
typical N- vortex configuration, with centers X^, (j=l, N), and windings rij is 

N 

(x - X ) 

0j(x,Xj) = arctan[-^ . (10) 

[x — y^jji 

The contribution of (|T0|) to the action is then easily evaluated 

vr ^ A 271 ^ 
Sn = — J2^j^k^og \Xj-Xk\ . (11) 

j=i ^ j<k 

In the last equation, the neutrality condition for the vortices' charges must also be 
assumed in order to have a non-zero action when the IR cutoff is removed. When 
also spin waves are present, the total action becomes the sum of (|Tl|) and the usual 
spin wave part. In what follows the (usual) approximation is made of considering 
only windings equal to ±1 for the vortices, since they are the most relevant. Due to 
this constraint, each sector must contain a fixed number of vortex-antivortex pairs, 
and we shall include a fugacity rj plus a combinatorial (classical) factor to take 
into account the indistinguishability of the vortices of equal charge. The partition 
function is then calculated by summing (with the appropriate weight factors) over 
all the possible configurations within each topological sector, and then over all the 
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topological sectors. Sumarising, 



oo 2N „ N 



^ = E 7^ / \{dX'^dYf Z^{{X,,Y,}) , 
iv=o ■) •' i=i 



(12) 



where Zj^ is the partition function for the spin variables in the presence of vortex- 
antivortex pairs, with coordinates X^, Yj-, respectively. One easily realises that the 
spin- wave contribution factors out, and the total partition function becomes the 
product of the partition function of a free scalar field by a Coulomb gas partition 
function. This Coulomb gas partition function can then be mapped to a Sine- 
Gordon partition function f^. However, the correspondence between the fields in 
both partition functions is somewhat ad hoc in this framework. 

2 The polynomial or 'first order' formulation 

We know from the previous section that the only non-trivial correlation functions are 
the 0(2)-invariant ones, and this is tantamount of invariance under translations in 0. 
It seems then natural to look for the possibility of describing the system completely 
in 0(2) invariant terms. We realize that the simplest possible ^-translation invariant 
field variable is 



where g is a constant with dimensions of mass, introduced to make dimensionless. 
We may of course rewrite the action (^) in terms of only, but we must also take 
into account the fact that is not entirely arbitrary, but a 'pure gauge' field, i.e.. 



When vortices are present, (|T^ is relaxed, allowing for a discrete set of points where 
the rhs is non-zero. 

We construct the first order action for the system in the spin wave sector by 
rewriting (^ in terms of and then adding a Lagrange multiplier term for the 



(13) 



ei^udf,L^{x) = . 



(14) 
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condition (|T^): 

(fx{—L^L^ - igOe^ydf^Lu) (15) 

where 6' is a scalar field which enforces condition (|T¥). 0(2) invariant correlation 
functions are then calculated in terms of L^, observing that 

exp{i[0(xi) - 0(x2)]} = exp[i5( / dz^L^{z)\ (16) 

Jx2 

where the line integral is taken along any smooth path joining X2 to xi. (|T6|) defines 
the Ihs in terms of the rhsQ. In order to give a consistent definition of a local spin 
field, we should require (|T6D to be path-independent. This will happen as long as: 



g f dz^L^iz) = 27tN (17) 

(where N is any integer) for every closed curve C. 

The general non-zero correlation function of the usual formulation will satisfy 
condition (H), then for each ej = +1 in (|^), there must be an = —1. Whence 
the general correlation function of the usual formulation can be constructed by 
forming 'neutral pairs' like the Ihs of ([T6|) , and then writing them in terms of the 
corresponding Wilson line on the rhs (of course there are many different ways to 
chose the pairings, all giving the same result). 

The generalization of ([T^ to the case when vortex-antivortex pairs are present 
is just 

e^yd^L^{x) = pn{x) 
2_ N 

Pn{x) = -Y.[S{x'X,)-5ix-Y,)], (18) 
9 j=i 

where we use the notation of Equation (12). Within the topological sector defined 
by (pISD , the generating functional of correlation functions is 

Zjv(J) = I VL^Ve exp(-5^ + j d\j^L^) , (19) 



"'^This is done just to make contact with the usual description, but there is no necessity to 
introduce the spin fields, neither to solve the constraints on within the polynomial formulation. 



where 

(fx {—L^{x)L^{x) - ig6 [e^^d^L^ - Pn{x)] } (20) 

with Pat as defined by (p^. Let us calculate then the spin-spin correlation function 
in this sector 

(s(x) ■ s{y))N = VL^Ve exp[z j^^ dy^L^{y)] exp(-5^) (21) 

(where 3? means real part of). Thus (21) can be obtained from Zn{J) just by 
specifying the current which reproduces the 'Wilson line' for L^, i.e., 

J,{u) = zg f'ds^6{u-z{s)) . (22) 
JO as 

We first calculate Zn{J) for arbitrary J and then we take it to be equal to (^2]). 
Integrating in (0), we get 

Zn{J) = jve eM-\d^.Od^O - ^ J2 + e{--^e^,d^J, + igpN)] ■ (23) 

The integral over 9 is also Gaussian. The final result is 

ZnU) = exp[ — ^ / (fxd ■ Jd~^d ■J+—f ^xpnO'^pn 
2g^ J 2t J 

- J d^xe^,d^J,d-^p] . (24) 

Using for the explicit form (^2]), ( p^ yields 

t vr ^ 

{s{x)-s{y))N = exp[-— log I x-y |]exp[-- ^ log I Xi-y^ I] 

AT 

X cos{^[a(x - Xk) - a{x - Yk) - a{y - Xk) + a{y - Yfc)]} 

k=l 

= G{x,y;{Xk,Yk};t) (25) 

where we have absorbed the divergent factors in a renormalization of the fields, and 
we use the function G to denote explicitly the dependence of the correlation function 
also on the coordinates of the vortices and the temperature. 
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Now we particularize equation (25) for the case = 1 

t 2tt 
G{x,y;X,Y;t) = exp[-— log | x - ?/ |] exp[ — - log | X - F |] 

ZTT t 

X cos[a{x-X)-a{x-Y)-a{y-X) + a{y-Y)], (26) 

where a{x) = arctanx2/xi for any x = {xi,X2). The correlation function (P^) is 
symmetric under the interchange of spin and vortex coordinates, plus a transforma- 
tion of the temperature: 

47r2 

x^X,y^Y ; t^ — . (27) 

Of course, this symmetry is inherent to the model, and not a consequence of using 
the polynomial formulation. However, we will need it to realise this symmetry as the 
result of some invariance of the path integral under a transformation of the fields. 
To achieve this, we need a description where spins field and vortices field appear in 
a more symmetrical way. is the spin field and one can realize that 6 plays the 
role of a vortex field. Indeed, the path integral representation of (p6|) as defined by 
(21) is 

(s(x) ■ s(y)) = 3^1 VL^Ve exp[ J £xi-^^L^L^ + igOe^.d^L,)] 



X eM^g r dz^L^) exp[-27r2(^^(X) - e{Y))] , (2^ 



where one can see that the exponential of 9 creates vortices when averaged with the 
spin wave action. We note that 



r'0exp[ j (fxigOe^ud^L^] exp[-27ri(0(X) - e{Y))] 

X exp[-2mg dz^9^{z)] , (29) 

where A is a new Lagrange multiplier field, which enforces the condition e^j^yd^dy = 
0, which is solved by 6^{x) = ^d^9{x). Inserting ( [29| ) into ( pH]) we get a more 



symmetrical description, in terms of the spin field and vortex field 
{s{x)-s{y)) = VL^Ve^VA 

X exp{ig f dz^Lf,{z)) exp[-2nig f dz^9^{z)] . (30) 

Jy JY 

Thus, the spin-spin correlation function corresponds to the average of the product 
of two 'Wilson lines', one for and the other for 6^: 

rx rX 

(s(x) ■ s{y)) = ^{exp[ig / dz ■ L] exp[2'iTig / dz ■ 6]) , (31) 

Jy Jy 

where the average is performed with the action 

2 

S = I d'^xi^L'' + igh^^e^K - tgAe^^dM ■ (32) 

Then, the duality transformation amounts to performing the following transforma- 
tion in thew action (IST 



1 47r2 
L^{x) ^ 27r9^{x) , 9^{x) ^ — L^(x) ,t^ — . (33) 

Zn t 

One easily verifies that the average of the same Wilson's lines operators with the 



transformed action gives the transformation (|27D . Note that the effect of this change 
of variables is to interchange the roles of vortices and spins, as well as low and high 
temperatures. In this sense, t = ^ is a kind of self-dual point, where spins and 
vortices are interchangeable. 

3 Mapping to the Sine-Gordon Model. 

The total partition function is defined following (|12|) 

oo 1 „ N 

^ = E 7T^2V''' / n d'X,d%Z^{0) , (34) 
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with 



2 

Z^(0)({X„rj) = I VL^Ve exp{- j d''x[^^L^L,-tgeie,,d^L,-pM)]} , 

PA^(a^) = - E ['^(^ - - '^(^ - • (35) 
9 j=i 



Performing the Gaussian integration over in (^5]), and using the delta func- 
tions in the definition of pn, we can rewrite Zn as 

N 

Zn = iU exp[2m{9iX,) - 9iY,))] ) , (36) 
i=i 

where the average () is defined by 

{F{e)) = JveF{e)exp[-s{e)], 

3(9) = ^- j (fxd,ed,e . (37) 



Using ( p7D in (|3^), we can rewrite the total partition function Z as 

^ = £ [/ ci^^e-^-^(^')]''[ / d^e^^-<^^y^{ ) . (38) 

Taking now into account the fact that only products of exponentials that satisfy the 
'neutrality' condition are non-zero, it is straightforward to check that 

Z = (exp(2r/ J d^x cos{2ne))) (39) 



or 



Z = jve exp{- J d^x[^{def - 2r/cos(27r^)]} (40) 

which is the desired Sine-Gordon action in terms of the Lagrange multiplier field 6. 
This coincides with the result obtained by more traditional methods. 

The effective action for 6 in (^) was obtained under the assumption that the 
vortices can only have charges equal to ±1. Let us consider now an extension of the 
model. It consists of relaxing this constraint, allowing for the charges to be equal 
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to plus or minus any real number q. Of course, the local spin interpretation will no 
longer be true, since the definition (^) requires condition (p!7|). In this sense this can 
be considered as an extension of the 0(2) model, which allows for vortices of non- 
integer charge. However, we will first assume a unique value for q in the model. The 
topological sectors will then have 2N vortex-antivortex pairs as before. Again only 
neutral combinations will have a finite weight. Of course we will also average over 
the positions of the vortices, and include the corresponding combinatorial factors. 
One can verify that instead of obtaining the Sine-Gordon partition function (^0]), 
we get 



Thus, even when the change performed on the spin system is drastic, the Sine- 
Gordon parameters change smoothly from g = 1 to any q. Things change more 
dramatically if we average now over all the possible values of g, since this produces 
delta- functions of the field 6. The result is a kind of field-theoretic delta-function 
model: 



4 Strings of vortices 

We have shown how the usual point-like singularities (vortices) are introduced in the 
polynomial formulation. Let us consider now string-like singularities (which could 
be regarded as strings of vortices). The obvious generalisation of the procedure we 
followed for the point-like case is to impose on a constraint like: 




(41) 




(42) 



(43) 



where pst is the density appropriate to a string: 




(44) 
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where the string's path is parametrized by 7(5) : [0, 1] —>■ R^, and q{s) measures 
the density of vorticity along the curve. Of course we can consider more than one 
string, just by using in the rhs of ( ^41) the sum of the densities corresponding to the 
paths. One can see that the spin configuration which corresponds to a state like the 
one defined by ( ^31) and (|4^) is: 



s[x) = [ cos (p[x ), sm (p[xj j , 

(j){x) = / (isg(s) arctan(x — 7(3)) . (45) 
Jo 

As for the vortices, the non-local definition of the spin field will be consistent only 
if condition (p^ is met for any curveQ. When strings are present this implies: 

fdsq,{s) = N, , Vj, (46) 

•J 

where the A'j's are integers. Also, in order for the action to be non-zero when the 
IR cut-off is removed, a neutrality condition must be satisfied 

TV 1 

E / ^As) = . (47) 

When only one string is present (^) supersedes (|i6|), and thus there is not true 
singularity since ( ^5|) implies that the net rotation of the spin in the local frame is 
null. Note that there is a crucial difference between open and closed strings, because 
while in the former the Wilson line definition works globally, in the latter it applies 
only to one of the simply-connected regions into which the space becomes divided 
(we recall that we cannot cross a singularity with the Wilson line). Let us estimate 
the Boltzmann weight of a configuration of a string-antistring configuration. It 
is straightforward to calculate the action due to this configuration in the polyno- 
mial formulation. It becomes just the Coulomb energy of the corresponding charge 
distribution: 

N 



Sst = - / ^-^1 / '^'^2 qi{si)qj{s2) log I ai{si) - l3j{s2] 



(48) 



^We do not allow the curves to intercept any singularity. 
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where a and parametrize the strings and antistrings paths, respectively (we have 
not written the (divergent) self-energies). With the natural definition of the mean 
position (Xj) of the j-th string, 

log(x - X^) - ^ dsq^{s) log(x - a{s)) , (^ dsq^{s) = +1) , (49) 

and analogously for an antistring, the action looks exactly like the one of N vortex- 
antivortex pairs, i.e., 

N 

Sst^ -Y. log \X3-Yk\ . (50) 

Thus we have arrived at the conclusion that the Botzmann weight of these config- 
urations is like the one of the usual point-like vortices. However, large strings are 
suppressed in the partition function because of the restriction about their positions. 
For example, the volume inside a closed string cannot be occupied by another one, 
and so they should be more strongly suppressed than the open ones. 
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